Exploring mudbrick architecture and its re-use in Artaxata, Armenia, during the 1st millennium BC. A multidisciplinary study of earthen architecture in the Armenian Highlands

Mudbrick constructions are extremely common in ancient western Asia, including the 1st millennium structures of the southern Caucasus and Armenian highlands. However, in the Caucasus the geoarchaeological study of these materials to provide insight into building practices and social structure is a topic little researched, especially when focusing on the longue durée. Artashat/Artaxata (Ararat region, Armenia) was the capital of the Armenian Kingdom of the Artaxiads, founded in the eighties of the 2nd century BC, but even before this the site was occupied in the Chalcolithic period, (ca. 5200–3500 BC), Early Iron Age (ca. 1200–900 BC) and in the Urartian period (ca. 800–600 BC) as well. All the previous occupation phases showed communities that made extensive use of earthen constructions as determined during past and recent archaeological excavations. This multidisciplinary study seeks to examine mudbrick architecture as a proxy for environmental and social interactions during the 1st millennium BC combining geoarchaeology, archaeobotany and building archaeology. We analyzed changes and continuities in architectural form and practices, alongside reconstruction of technological and social processes, to identify issues of raw material procurement, attestation of re-use, and consistency of building practices. The results of the geoarchaeological analysis of the earthen building materials used in different parts of the ancient city point to a re-use of materials over time.


Introduction
Compared to Mesopotamia and the Levant, mudbrick constructions in the Armenian highlands and Caucasus have been the focus of limited research, mostly at the macroscopic and chronological level [1][2][3]. The importance of architecture, especially earthen architecture, in shedding light on past social practices and interactions is well attested [4][5][6], both as a proxy to understand past environmental changes and to determine past social practices [7,8]. Nevertheless, Classical period buildings are often not given the same attention as prehistoric ones in the analysis of earthen building materials and their labor organization. The Armenian-German Artaxata Project (AGAP) launched in 2018 at Artaxata/Artashat (Ararat region, Armenia) focuses on the urbanism of the capital of the Armenian Kingdom of the Artaxiads, founded in the eighties of the 2nd century BC. In Armenia, specifically Artaxata, the Classical period ranges from the 2nd century BC to the 4th century AD. This project builds on the results of the previous Armenian excavations of the 1970s [9][10][11], which have studied the urban development of the numerous hills that constitute the ancient city (Fig 1).
Our research first aims at bringing forward the analysis of earthen building materials at the site of Artaxata, located in the Ararat Valley of central Armenia. Second, we compare construction practices between the two main phases of site occupation during the 1st millennium BC: Urartian and Hellenistic. This comparison allows us to investigate differences and similarities in sources of raw material procurement and building practices, but especially to provide a new methodology to assess the issue of re-use of mudbricks from archaeological context. This is a particularly important question, as currently there is no such methodological framework for the determination of earthen building material (henceforth EBM) re-use.
The re-use of EBM, specifically mudbrick, is an essential and understudied topic in archaeological research. Their re-employment in buildings over the centuries does not follow the same pattern of other spolia and ceramic building materials such as laterizio. Re-used mudbricks are taken from previous buildings, broken down and remixed with additional water and temper to mold new bricks. In contrast to fired bricks in which re-use is easily detectable, only rarely are mudbricks taken from a building to be directly employed in another construction (author personal observation during her work in Egypt). More commonly they are worked and reshaped thus fully recycling the materials, but also re-initiating the manufacturing process [12]. This type of intensive recycling can still be considered re-use with the caveat that it is incredibly challenging to trace and identify in the archaeological record.
The traditional mudbrick chaîne opératoire consists of collection of raw materials, thus soil, water, and human-induced inclusions used as degreasers such as the more commonly used vegetal temper, but also shell fragments, sand, and crushed plaster. The mixing, molding, and drying activities have been extensively described [13,14], indicating how the analysis of the chaîne opératoire can provide useful information regarding the procurement of raw materials within a community [8,15], the presence of skilled or unskilled workforces [6,16,17], and the development of apprenticeship within a community [14]. The follow up study of how earthen structures decay, collapse and participate in postdepositional processes brings forth further knowledge regarding the cycle of architecture and materials collapse [18,19]. Here, we argue that an integrated approach that combines geoarchaeology, micromorphology and archaeobotany provides us with the necessary parameters to recognize re-used building materials.
One innovative aspect of the present study is the phytolith analysis of EBM thin sections. While the observation of opal of biological origin within petrographic thin sections is not uncommon [20,21], the systematic examination of phytoliths in ceramic thin sections is restricted to a few case studies [see for instance [22][23][24][25][26][27]. Even more infrequent are studies elaborating on the phytolith content of construction remains such as bricks and cement [18,[28][29][30][31][32][33]. Most of these studies extract phytoliths which ends in the loss of part of the phytolith's context.
Phytolith taphonomy is a complex matter [for a discussion of this issue see [34][35][36]. Prior to any taphonomic process such as erosion, corrosion or transport, decomposition of organic matter needs to occur to release phytoliths. Decomposition of the organic matter can occur either before or after the incorporation of the plant fragment into an archaeological record. When decay of the plant fragment takes place after its incorporation into the archaeological object, the relative distribution of the phytoliths, as observed while they are in anatomical position, should be preserved, provided that no post-depositional perturbation occurred. This relative distribution, together with its context, is destroyed when phytoliths are removed for analysis using typical extraction processes. On the contrary, they are preserved in thin sections.
Phytolith analysis on thin sections is a technique that inventories the distribution patterns of the phytoliths within archaeological contexts as well as the phytoliths composing each of these patterns. Each pattern (namely, Isolated, Clustered and Articulated (for a definition see [37,38]) hides a different history. The inventory of each distribution pattern together with the phytolith composing each of them, enables to discriminate phytoliths sharing (or not) a common (post)depositional history and botanical origin [37,39,40]. Concerning ceramic thin sections, it discriminates phytoliths observed in voids from those in the clay matrix. Previous taphonomic research determined phytoliths within voids derive from plant material selected as temper while those in the matrix relate to the natural composition of the clay [37]. Such an approach successfully contributed to the identification of ware source areas and the reconstruction of the chaîne opératoire [41].

Archaeological background
The Hellenistic city Artaxata was founded in the 180s BC by the Armenian king Artaxias I as the capital of his kingdom at a location which was said to be previously uninhabited [42,43]. However, recent archaeological research has shown that the city was the site of a considerable Urartian settlement and Artaxias I resettled the apparently abandoned site [44].
The Hellenistic city was located on 17 natural limestone hills which rise from the Ararat Valley and its lower city stretched into the plain. The river Araxes passes the west and south of the city. The rivers Metsamor and Azat, which flow into the Araxes, come from the northeast, while the Vedi waters the south of the city. The Ararat Valley is a rich and fertile alluvial plain with clay deposits used for pottery and mudbrick production [45]. The capital city of Artaxata soon became a major metropolis in Armenia and a melting pot of a variety of interactions with the Mediterranean, the Levant, Iran, Mesopotamia and the northern Caucasus. Architecturally, Artaxata was strongly reliant on local construction techniques and most of the buildings were made of mudbricks placed on low stone foundations. In the archaeological record, at least the lower courses of the mudbrick walls can be found intact, even if most of the former walls are eroded and deformed.
Armenia was affected heavily by the antagonism between Rome and Parthia but especially under Tigran II was a mighty transregional power. The city experienced several violent destructions, among them around 66 BC a partial destruction during Tigran the Younger's rising and in AD 59 by Corbulo. These destructions seem to be attested in the archaeological record. In this study mudbrick samples are included which stem from the Urartian period and from the first two main phases of the Hellenistic settlement, from Phase I (180s to approx. 66 BC) and Phase II/III (approx. 66 BC to AD 59). These mains phases were identified by stratigraphic excavations and a 14C sampling program. The data was mainly assembled on Hill XIII which is the main site of the Armenian-German Artaxata Project (AGAP) (Fig 2).
Although our research focuses on the Hellenistic Lower City and its three-phase development (Phases I, II and III) during the Hellenistic period, Artaxata's nature as a multi-period site is pertinent to our study. Besides the dense Hellenistic settlement, traces from the Chalcolithic, Middle Bronze Age, and-as previously mentioned-Urartian period have been documented [46]. During this work, remains of an apparently extensive Urartian settlement were found on Hill XIII and in the southern plain, alongside the previously documented citadel on Hill II [9][10][11].
More importantly the radiocarbon dates provide us with solid evidence that these two distinct periods of occupation, Urartian and Hellenistic, occurred in the same area of the settlements indicating for instance that Trench 11 areas present Hellenistic occupation built directly on top of Urartian remains, even if the exact stratigraphic sequence needs to be further investigated (Table 1. Data in S1 File).
Geological background. Armenian geological map is complex and multifaceted. Most of the area around Artaxata is characterized by diatomite clays and Quaternary lacustrine alluvium (Fig 3). Paleogene and sandstone deposits alternating marls, sandstones, siltstones, tuffs and gravel surrounded the site to the North and the North-East. Directly to the East of the site, the geological landscape includes cretaceous deposits of metamorphic schistose limestones, tuff conglomerates and basalts, while in the south we have cretaceous carbonate and limestones. Ting et al. conducted a sampling of the area to study the raw materials for ceramic production in Artaxata, which is the baseline for the geological analysis in this article [45].

Excavation area
The main focus of the AGAP, starting in 2018, has been on Hill XIII and the adjacent plain to the south. The hill is located at the transition from the Upper to the Lower City of Hellenistic Artaxata. The starting point for the excavation work is the results of geophysical prospection, which revealed several building structures in the Lower City [47]. In five excavation campaigns, Hill XIII, an anomaly in the form of a dotted line to the north of it, and a hall-like structure recognizable in the magnetics in the southern plain were investigated.
While the anomaly to the north turned out to be the remains of pillars of an unfinished Roman aqueduct bridge [48], the remains of large-scale constructions were found on Hill XIII. The most significant features, which were to a certain extent previously evident in the magnetics, are quarry stone wall bases, some of which overlap, revealing a relative chronological sequence. Almost the entire eastern summit of the saddle-shaped Hill XIII and its slopes as well as the depression on the western hilltop are covered with a massive layer of collapsed mudbrick, deriving from the walls that were built on the quarry stone bases [49].
Construction here took place in two main phases, which could be firmly dated using 14C samples. Phase I begins around 180 BC with the construction of a presumably sacred building with a broad central room, furnished with rich stucco decoration [48]. Traces of fire and partial destruction were recognized, which occurred around the middle of the 1st century BC, not impacting the mudbrick walls. In the sub-phase of Phase I (Ib), a change in function to domestic architecture becomes apparent. This change in function becomes clear in Phase II (from the middle of the 1st century BC). So far, two corridor houses separated by an alley have been excavated from this phase. Phase II was followed by the almost congruent Phase III. Structural and functional changes can hardly be observed between the two phases. The end of Phase III is marked by the destruction under Corbulo in 59 AD. The small number of finds allows the assumption that the inhabitants had already left their dwellings before the destruction. Hill XIII was then never resettled but visited only sporadically and the dwellings left behind weathered away. Already in 2018, remains of past mudbrick architecture were documented in a test trench in the area adjoining the eastern hilltop to the north [50]. When work was resumed in 2022 it turned out that we are dealing with a massive filling underneath a retaining wall of Phase I. According to the pottery finds the material used for leveling this area must come from the Urartian period.
Already in 2018, a structure with a length of over 60 m was discovered on the magnetogram in the plain in the southern part of the investigation area [44]. This was partially excavated in 2021. In an area of 195 m 2 the northwest corner of the structure and a three-part double gate in the center were uncovered. The function of this structure is still unclear, but on the basis of the 14C data, it can be dated with certainty to the Urartian period around 800 BC. This represents an important finding in pre-Hellenistic urban history and the first monumental complex from that period in the study area.
In addition to the abovementioned features and some scattered finds, a massive mudbrick wall on the northern slope of Hill II is a particular testimony to the Urartian past of the site. It is assumed that the remains still visible today were part of an Urartian fortress with a size of 2.5 to 2.6 ha [51] (Fig 2).

Materials and methods
Sampling took place from different sectors of Artaxata that have been recently excavated and from the Urartian defensive wall, excavated in the 1970s. When possible, the samples were selected from multiple rows, top to bottom, of the Urartian defensive wall (Hill II), preferably from external walls of the buildings identified on Hill XIII, Trench 11 in the plain, which presented both Urartian and Hellenistic levels, and trench 1-4, 6, 8, which displayed Hellenistic occupation levels (Fig 4). Due to the continuous rebuilding and levelling actions during the Hellenistic period, often only limited mudbrick remains such as one row of mudbrick, were preserved in situ for sampling. Therefore, the 34 samples include the multiple phases of site occupation, specifically Urartian and Hellenistic, including phase I, II and III registered on Hill XIII ( Table 2).
The 34 samples underwent geochemical and petrographic analysis. Based on the results and material preservation, we sub-selected 17 samples to perform granulometric analysis, calcimetry, loss on ignition and four for phytolith analysis.

Methodology
Major, minor, and trace element concentrations of 34 bulk samples were analyzed with a Rigaku NEX-DE VS bench-top ED-XRF spectrometer housed in the University of Helsinki Laboratory of Archaeology. The instrument was operated in point analysis mode, using 1 mm beam diameter and a camera view to select inclusion-free points to acquire paste geochemical compositions. The results are normalized means of 3-5 analyzed points per sample, measured in helium atmosphere, using a tube voltage of 60 kV, 35 kV, 6.5 kV and acquisition times of 60, 60 and 100 s for high-Z, mid-Z and low-Z elements, respectively. The results were quantified by the instrument's software and fundamental parameters; standard reference materials Burnt Refractory NIST 76a and Brick clay NIST 679 were analyzed with the samples to control data precision and accuracy. The ED-XRF dataset was processed with the IBM SPSS 28 software, the compositional groups are based on cluster analysis (CA) groupings using the Centroid Clustering method run with the elemental concentrations of MgO, Al 2 O 3 , SiO 2 , K 2 O, CaO, TiO 2 , MnO, Fe 2 O 3 , NiO, ZnO, Rb 2 O, SrO, Y 2 O 3 , and ZrO 2 detected in the entire sample set (P 2 O 5 , SO 3 , and Cl were excluded due to suspected contamination).

PLOS ONE
Exploring mudbrick architecture and its reuse in Artaxata, Armenia, during the 1st millennium BC

PLOS ONE
Loss on Ignition (LOI) was used to quantify organic matter content as the samples were pre-dried for 24 hrs, cooled, and weighed. Then samples were placed in the furnace at 500˚C for at least 6 hrs, cooled and weighed again for organic loss.
CaCO 3 was calculated through an automated calcimeter, GEO-RS calcimeter, capable of measuring the percentage of carbonate from the pressure increase within a sealed circuit based on the carbon dioxide developed from the reaction of 0.500 g of sample with 5 ml of 10% (v/v) HCl. The calcimeter was calibrated using 0.500 g of 99.95% pure calcium carbonate (produced by Alfa Aesar) and 5 ml of 10% HCl. The calibration was repeated after every 5 analyses performed.
Petrographic and mineralogical analyses of the thin sections were conducted with a Leica DM2000 polarized light microscope with an attached digital camera, working with a magnification between ×5 and ×40. The samples were analyzed and grouped following the system of structure and component descriptions [for detailed methodology see [52][53][54][55][56].
A few samples were selected for micromorphological analysis. They were oven dried at 50˚C and impregnated with polyester resin without disturbing their original structure. When they were solidified, they were cut with a rock saw into slabs of 2 cm thickness and with gradual cutting and polishing, thin sections of 30 μm thickness (5×7cm) were created. The initial processing of the samples was conducted at the M.H. Wiener Laboratory of Archaeological Science, ASCSA, while the thin sections were produced at Quality Thin Sections (Tucson, Arizona). In total, four thin sections were produced. The thin sections were photographed in high resolution [see 57]; they were studied under the stereoscope and the polarized microscope in magnifications ranging from ×1.25 to ×40 [55,58].
In addition, four petrographic thin sections were also used to study phytoliths and testing the viability of the methods for earthen architecture. Observations were conducted on a Zeiss Aksioscope under PPL and XPL [59]. When needed, additional observations were conducted under UV and Blue lights. Thin sections were scanned along two horizontal and vertical lines at magnifications ×100, ×160, ×200, ×500, and ×800. The purpose of the scanning at low magnifications (×100, ×160, ×200) is to detect the presence of phytoliths in voids. Scanning at higher magnifications (×500 and ×800) aims either at naming phytoliths observed in voids or detecting phytoliths in the fired clay. Indeed, due to a potential wedging effect [58], phytoliths can be masked by the fine fraction. Naming of phytoliths follows ICPN 2.0 (ICPT [60,61]). As recommended for the analysis of soil thin sections [62], attention is also paid to the presence/ absence of opal residues of biological origin (diatoms, chrysophyceae cysts, sponge spicules, external amiboid skeletons).

ED-XRF results
The cluster analysis dendrogram (Fig 5 and Table 3) of the ED-XRF geochemical dataset indicates one main cluster (Cluster 1), including 26 of the 34 analyzed samples, and compositional outliers which fall outside the main group. Cluster 1 is a chronologically mixed group of calcareous mudbricks (CaO c. 12.5-20 wt%), characterized by relatively high silica values (ca. 49 wt% on average) and iron values at ca.10 wt% on average. Of the eight data outliers, samples AA1, AA7, AA12 and AA34 display elevated calcium, iron, nickel, copper, and zinc oxide values; the NiO concentration of AA1 at almost 3000 ppm implies nickel contamination. Furthermore, samples AA8, AA10, and AA22 all display exceptionally high sulphur concentrations (18-24 wt%), accompanied by cobalt oxide levels at c. 120-230 ppm, potentially linking these sample materials with volcanic activity. It is also noteworthy that sample AA21 clusters with the main group, yet also shows enriched SO 3 levels at c. 11 wt%. In addition, potential exposure to metals is suggested by the enriched nickel value of sample AA6 (NiO at c. 900 ppm), the relatively high copper content of AA17 copper (CuO c. 400 ppm), and the elevated zinc, copper and cobalt values of AA4. In addition, sample AA15 shows the highest Sr, Zr and Ba oxide concentrations measured in this sample set (Full ED-XRF results are provided in S2 File).

Petrographic and micromorphological results
The results of the petrographic analysis point towards one main fabric, matching the results from the geochemical analysis. This fabric is characterized by the presence of igneous rocks in a greenish-brown groundmass (Fabric 1, n = 30), with an outlier group formed by a lone individual (Fabric 2, n = 1).
Petrofabric 1 is very homogeneous in terms of geological composition, although the distribution parameters allow for a sub-classification of individuals (Sub-fabric 1.  rocks]) (Detailed petrographic data in the Table in S3 File). The grain size of the aplastic inclusions is generally bimodal, densely packed in the fine section and with a coarse fraction single to double spaced ranging from medium silt to granules (0.3-3 mm). The main characteristic of this petro-group is the presence of polycrystalline sub-angular to sub-rounded extrusive rocks (andesite and basalt) together with pumice fragments and greenish-yellow volcanic glass. Aplastic inclusions also include other types of materials, but are less frequently present, in a poorly sorted grain distribution with angular to sub-rounded sphericity of the particle. It is the case of calcite, plagioclase feldspar, mono and polycrystalline quartz, metamorphic rocks (schist, phyllites), serpentine, biotite, iron-oxide nodules, pyroxenes (augite, hornblende), and rare limestone. Almost all samples of the group present a few mudstone grains of variable texture and random size, including a particular type of grey mud with chert and brown silt loam lumps. The presence of microfossils is constant and varied but not frequent. We have identified Planktonic and Benthic foraminifera, Echinoids, ostracods and gastropods. The presence of micro and mesovoids is common, but not excessively frequent, and they mainly follow the shape of vesicles or channels. The relationship between these channels and the use of vegetal degreasers is clear in the individuals AA1, 2, 3, 4, 8, 17, 20, 25, 27, 30, 31 and 34. Last, it has also been possible to identify other materials as added temper, such as charcoal in samples AA2, AA4, AA7, AA10, AA17, AA18, AA20, and AA34, sometimes reaching a large size and preserving its microstructure.
Three thin sections were produced from sample AA34 (a-c) dated to the Urartian period and one from sample AA8 dated to the Hellenistic period. These samples are the biggest and most diagnostic samples collected in order to identify deformation features.
The AA34 samples include calcareous yellowish-brown silt loams, moderately sorted with subrounded to rounded and polyconcave voids, reaching 1 cm in diameter; elongated channels are also recorded, locally curved, reflecting plant imprints of maximum 1 cm size, in irregular distribution and orientation or spiral deformation (Figs 7, 8A and 8B) (Detailed micromorphological description in S4 File). The voids are locally banded, parallel oriented and distributed, indicating shear zones (Fig 7A). Porosity is high, reaching 20%. The b/f fabric is speckled to crystallitic or undifferentiated (Fig 8D). Sand grains are sub-rounded to sub-angular and include quartz, plagioclase feldspar, chert and igneous rocks. Pumice grains are also included (Fig 8C).
Subrounded dark yellowish-brown calcareous aggregates are randomly distributed (Fig 8B). In sample AA34c these aggregates are very pronounced (Figs 7B and 8). Dispersed microcharcoal and charcoal fragments, along with calcitic sediments suggest the presence of ashes (Fig 8D).
Moreover, aggregates with laminated microstructure indicate the addition of alluvial sediments (Fig 8F).
While these results matched the initial petrographic analysis, micromorphological analysis also indicates that both the calcareous aggregates and the matrix are related to the addition of lime in the sediments. Lime is identified microscopically by the presence of a dense calcitic cementing fabric with occasionally shrinkage cracks, a few vesicles, and irregular voids with smooth walls (Fig 8B). Deformation features in these samples include the aligned and parallel oriented vughs and voids as described above, which are produced by the compaction of wet sediments when air is trapped in the voids [63]. The curved orientation of fine plant imprints further shows the direction of the shear stress while also providing clear evidence of vegetal

PLOS ONE
tempering. The same pedofeature has been recognized in the Urartian Sample AA20 in thin section. Another indication of the same process is the rounded spiral deformation of plant remains, which is indicated by their imprints. In sample AA34b, deformation features are illdeveloped and limited to linear, wavy, and polyconcave voids. However, the most conspicuous feature is the presence of lime lumps in the form of poorly crystalline calcareous aggregated areas, with dark gray appearance and low birefringence, indicating only partial reaction and carbonation [64,65].
The Hellenistic sample AA8 is largely identical to AA34. It therefore includes yellowish brown calcareous silt loams (Table in S4 File), moderately sorted with subrounded to rounded voids and elongated channels randomly distributed, locally curved, indicating plant imprints in irregular distribution and orientation (Figs 7C and 8F). The b/f fabric is speckled to crystallitic or undifferentiated. Crystallitic b/f fabric is tentatively attributed to the presence of ashes and spherulites are visible, indicating that dung was likely added to the matrix as a possible degreaser or has accidentally been included in the raw materials.
Porosity is still high (20-30%). One important novelty is the fine-grained and calcareouslime aggregates identified in AA34, which are also recorded here although in lower proportion.

The results of granulometric and calcimetric analyses
The results of the granulometric analysis suggest a slight differentiation in recipes between the Urartian and Hellenistic mudbricks (Fig 9). Urartian mudbricks seem to have a more consistent silt fraction-on average more than 50%-while Hellenistic mudbricks are characterized by a fabric with a silt fraction equal to or smaller than 40%. The organic percentage calculated through LOI is quite consistent between samples ranging between 2.5% and 5%. On the other hand, the percentage of CaCO 3 equivalent shows variation between individuals of the Urartian and the Hellenistic period. Calcium Carbonate is higher on average in the mudbricks of the Urartian period, although its presence is also attested in Hellenistic mudbricks. This variation may be linked to the use of lime as a secondary human-induced temper alongside chaff in Artaxata mudbrick production (Table 4). The PCA combining granulometry, Organic %, CaCO 3 % and selected geochemical data (Al 2 O 3 , SiO 2 , P 2 O 5 , SO 3 , Cl, K 2 O, CaO, TiO 2 , MnO and Fe 2 O 3 wt%) provide us with two different clusterings, highlighted by the confidence ellipses created by coding a multivariate t- distribution. The two ellipses underline a general communality of the raw sources but slightly different recipes between the Urartian and Hellenistic periods (Fig 10). The main difference, as also confirmed by the micromorphological and petrographic analysis, is the ratio variation between the fine and coarse fractions (Fig 11).  The results of Phytolith analysis The petrographic thin sections AA7, AA17, AA27, and AA34 were selected for exploratory analysis of the phytolith content. They represent a combination of Urartian and Hellenistic mudbricks from the two main occupation periods. Table 5 documents the absolute number of phytoliths captured by each of the observed distribution patterns in the fabric (Illustrations of the observed distribution patterns are provided in the tables in S5 File).

PLOS ONE
Exploring mudbrick architecture and its reuse in Artaxata, Armenia, during the 1st millennium BC By associating three distribution patterns (isolated, clustered and articulated) together with a more marked morphological diversity and the largest number of phytoliths, the fabric of AA7 differs from all the others. Basically, the AA7 fabric attests the richest and most diverse phytolith content. A silica skeleton was also recorded within it (Fig 12A and 12B).
The other fabrics provide additional specific features. AA17 is characterized by an association of isolated and clustered phytoliths where each pattern captures an almost equivalent number of phytoliths (of 9 [Isolated] and 10 [Clustered] respectively). Only isolated phytoliths were observed in the AA34 fabric (Table 6).
Analyzing the petrofabrics, we noted the presence of numerous phytoliths and diatoms. While phytoliths were present in all the analyzed thin sections, diatoms were not observed for AA34 but were for AA7, AA17, and AA27 (Fig 12C and 12D).
Even considering the limited number of thin sections analyzed, the AA7 petrofabric appears to be the one richest in phytoliths, the total number of observed phytoliths being at least 2.5 times greater than any others.

PLOS ONE
Exploring mudbrick architecture and its reuse in Artaxata, Armenia, during the 1st millennium BC

Discussion
The combination of results provides us with a new perspective on human-environment interaction in Artaxata during the 1st millennium BC, specifically during the Hellenistic period. The chemical and petrographic data point to a local source for procurement of raw materials for earthen building materials. Soil resources can easily be identified within walking distance from the site, including the geological deposits located within a radius of a few km from the site (Fig 3). For instance, the alluvial deposit that constitutes the backbone of Artaxata geology is the main source of soil procurement strongly influenced by the volcanic materials north of the settlement and from nearby Mount Ararat, which are characteristic of the matrix inclusions and aggregates in all of the mudbricks manufactured at Artaxata. The presence of polycrystalline extrusive rocks such as andesite and basalt in petrofabric 1 together with pumice fragments and greenish-yellow volcanic glass perfectly matches up with the geological

PLOS ONE
Exploring mudbrick architecture and its reuse in Artaxata, Armenia, during the 1st millennium BC characteristics of Artaxata's environment, very conditioned by the igneous materials from Mount Ararat [66,67].
On the other hand, the ratio of clay, silt and sand, carbonates and the phytolith inventory provide indications of multiple manufacturing phases and likely material re-use.
The PCA mapping of the granulometric analysis, organic percentage, chemical compounds and CaCO 3 equivalent highlights the presence of two distinct clusters (Fig 10). When we incorporate the analysis of the grain size and the matrix typology in thin sections, we noted a change of the mudbrick recipes between the Urartian and Hellenistic periods. Although some Hellenistic samples have been classified in the coarser matrix Sub-Fabric 1.1., a thinning of the recipe, which means a smaller fine fraction, characterized a significant part of the Hellenistic mudbricks and could be clearly associated with cyclical re-use of earthen building materials found at the site by the Hellenistic communities that worked on mudbrick manufacturing. The collection of new raw materials is often a time-consuming and economic burden for manufacturers [68]; thus, the presence of abandoned structures at the site could provide the necessary material while removing a consistent chunk of the traditional chaîne opératoire [14]. To test this hypothesis, we created multiple micromorphological sections to analyze the morphological differences between mudbricks of different periods.
The micromorphological studies support our initial hypothesis, identifying the aggregates in the Urartian sample AA34 as immiscible, suggesting the incomplete mixture of diverse materials (Fig 8B). In a few cases, sediments are characterized by a laminated microstructure, indicating their alluvial origin ( Fig 8E) and matching the geological deposits around Artaxata. More often the Urartian sample AA34 includes fine-grained and calcareous-lime aggregates, as described above.
This observation is tentatively interpreted as correlating well with the results of granulometry and LOI. More specifically, the lower percentages of silts and CaCO 3 in AA8 possibly reflect the lower number of fine-grained lime inclusions in comparison to AA34. The deformation features are here limited and ill developed with abundant polyconcave vughs. The matrix is crystallitic due to the presence of ashes and spherulites are visible, indicating that dung was likely added to the matrix as a possible degreaser or has accidentally been included in the raw materials. The presence of pedomorphic voids and the constant organic percentage indicates that organic matter was re-added during the Hellenistic manufacture and mixed with water and the broken-down mudbricks, a key source of appropriate mudbrick sediments.
In line with this argument, the presence of immiscible aggregates of sediments in Hellenistic sample AA8, may be associated with fragments of re-used broken-down mudbricks (Fig 7C), which presented characteristics clearly in line with the petrography of Urartian mudbricks (i.e., Vegetal temper and mineral inclusions), and stands out in the Hellenistic mudbrick matrix as not completely integrated into the new structure.
The data presented and the difference in morphological and grain size percentage between Urartian and Hellenistic samples seems to point to a specific process of re-use of building materials by Hellenistic communities. Urartian mudbricks standing in situ had progressively deteriorated due to abandonment, wind erosion, water erosion, temperature related deterioration, and chemical related deterioration [69,70], leading to a loss of the fine fraction in the bricks. While wind erosion tends to affect the coarse fraction within a mudbrick more, the other three types of deterioration affect the fine fraction equally. For instance, AA9 shows anomalously high sodium and chloride oxide levels at c. 20 wt% and c. 13 wt%, respectively, indicating possible salty water exposure or sodium chlorite mineral content. At Artaxata, the high values of soluble substances such as Na + , S, Cland Ca in the soil combined with the high porosity of the mudbricks and the extensive presence of water that permeated the bricks over time, such as in the case of precipitation, generated both water and chemical related erosion, thus deeply affecting the structural integrity and the fine fraction of the bricks [18].
A similar phenomenon is observed in many earthen structures when left exposed without continuous maintenance [18,70,71]. Mudbrick left exposed undergoes rapid deterioration that especially affects the fine fraction due to the loss of material [19,31], but also the cementing of smaller grains by salt crystallization such as calcium carbonate [71].
Once these were broken down and re-molded during the Hellenistic period, there was little to no addition of new sediments, mainly supplementing the recipe with new vegetal temper and other aggregates. This explains the same approximate percentage of organic temper and CaCO 3 between Hellenistic and Urartian mudbricks, but phytoliths clearly visible within voids were recorded only for AA17. This implies that plant material was used as temper. Their inventory tends to support the use of diverse plant material alongside by-products of cereal processing. The taxonomical occurrences of BLOCKIES are diverse and very common in the leaves of Cyperaceae and Poaceae as well as in dicots and conifers [61]. Hence, one cannot exclude the use of diverse types of plant material as temper, among which is the by-product of cereal processing.
Urartian mudbricks tend to show a consistently denser structure with a few thin elongated features where fibrous vegetal temper was present, as well as rotatory features as a result of the extensive kneading, On the other hand, Hellenistic mudbricks are characterized by a less dense structure with crudely aligned and rotated voids as evidence of limited tactile manipulation [14], in which the incorporation of previous EBM may be detectable.
Later processes could determine the phytolith content (distribution patterns, phytolith amounts and morphotypes) of the Hellenistic mudbricks. The phytolith inventory reveals that almost all phytoliths distributed along the Isolated and Clustered patterns are clustered within the Hellenistic petrofabric (Table 5). Clusters consist of a group of disarticulated phytoliths where not all phytoliths are necessarily of the same type or share the same orientation (Figs 12  and 13). These are understood as typical markers of post-depositional perturbations resulting from the physical reworking of the original material [59]. As such, the clusters observed for the Hellenistic mudbricks, as the studied material allows us to characterize them, could result from the re-use of Urartian mudbricks.
Granulometry, specifically the measurement of the ratio of clay-silt and sand, remains one of the most reliable variables in the analysis and classification of mudbrick recipes [4,6,12], but both micromorphology and phytolith analysis reveal themselves as key methodological additions to enhance our ability to discern cases of re-use. The new methodology employed, specifically the phytholiths' analysis and micromorphology, allows us to better understand the relationship between the natural and built environment, raw source-material procurement choices in the Hellenistic period, and the intensity of the construction process.
The recycling of mudbrick in architecture is not a new phenomenon [4,72], but in most cases the argument for re-use discusses the removal of finished artifacts, in this case bricks, from ancient walls. In our argument we discussed a different phenomenon in which previous building materials are re-used as an available and ready-made source of soil. This shortcut alters the traditional relationship between communities and the environment, indicating the opportunistic nature of Artaxata Hellenistic communities regarding mudbrick manufacturing and the chaîne opératoire. This may have been prompted by the need for a fast construction process and the necessity to build a new center in a short period of time.
Additionally, the benefit of using already existent earthen building materials stresses the rational use of sources of raw materials within an almost carbon-neutral construction process that cut the energy inputs, such as carrying the raw material up a steep hill, to maximize efficiency.
The evidence supports the written sources that, at the time of Artaxias I's foundation in the 180s BC, the previous settlement had been abandoned for quite some time and the previous structures had deteriorated. Since Artaxata was not continuously inhabited during the 1st millennium BC this seems to indicate that, during Hellenistic Phase I and II/III, the community re-used abandoned and exposed EBM in order to have economical and advantageous access to soil sources for their constructions.

Conclusions
Our research built upon previous mudbrick research, emphasizing the relevance of geochemistry and grain size analysis in the study of mudbrick, but also provide new evidence on the methodological importance of micromorphology and phytolith analysis to discern cases of earthen building material re-use. The data clearly indicated that an interdisciplinary approach is better suited to investigate the complex relationship between the architecture, natural environment, and past communities.
Finally, the recycling of earthen building materials is a well attested phenomenon that is often difficult to attest in the archaeological record. Our research presents an innovative methodology to observe reusing of mudbricks in archaeological sites employing geochemistry, micromorphology and phytoliths analysis. The method is relevant for archaeological and architectural studies, specifically those dealing with mudbrick architecture.